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Abstract. In addition to being suitable for laser cooling and trapping in a 
magneto-optical trap (MOT) using a relatively broad (~5 MHz) transition, the 
molecule YO possesses a narrow-line transition. This forbidden transition between 
the X 2 £ and A' 2 A 3/2 states has linewidth ~27rxl60 kHz. After cooling in a MOT 
on the X 2 T, to A 2 H 1/2 (orange) transition, the narrow (red) transition can be 
used to further cool the sample, requiring only minimal additions to the first stage 
system. The narrow line cooling stage will bring the temperature from ~1 mK 
to ~10 fiK, significantly advancing the frontier on direct cooling achievable for 
molecules. 
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1. Introduction 

Cooling of molecules has wide-ranging applications, from tests of fundamental 
symmetries mm, study of long range anisotropic interactions [1] , to cold chemistry 
01- The more complex structure of molecules as compared to atoms, though it 
provides a rich field of study, also makes difficult the cooling methods readily available 
to atoms. Indirect methods have so far produced the coldest molecules, via magneto¬ 
association and adiabatic transfer |7]. Unfortunately, the number of species that can 
be produced in this way is limited, and so we must look to direct cooling methods to 
support a wider range of molecular species. Cold molecular beams can be produced 
via supersonic expansion [8] or from a cryogenic buffer gas beam source [9]. Such 
beams can be decelerated with techniques such as Stark deceleration mm, Zeeman 
deceleration m, or via centrifugal deceleration m■ From a low velocity beam, 
a two photon process has been used to directly load a magnetic trap M- Once 
trapped, molecular collisions and evaporative cooling have also been demonstrated for 
the hydroxyl radical m- 

For a long time, it was thought that laser cooling of molecules was impractical 
due to the large number of states and therefore lasers required to implement an optical 
cycling transition capable of many photon scatters. However, in 2004 it was proposed 
m that a large class of molecules has sufficently diagonal Franck-Condon factors to 
be amenable to laser cooling. In 2008 the first concrete proposal was made for a class 
of cycling transitions that through angular selection rules avoids all dark states for 
molecules, enabling quasi-closed cycling and so the scattering of many photons [TT j. In 
the intervening time, optical cycling has been demonstrated in a number of diatomic 


Prospects for a Narrow Line MOT in YO 


2 


species, leading to demonstration of laser cooling (T 8 ], radiation pressure slowing 
m nn mi, and magneto-optical trapping in both two ( 2 2 ] and three dimensions 
[25], 1 25], Additionally, it has been proposed that the optical biclrromatic force [25] 
or ultrafast stimulated slowing [25] can be used for molecular beam deceleration and 
cooling m- 

The lowest temperature realized with these direct cooling techniques is so far that 
of the magneto-optical trap (MOT) [23], with temperature ~2 mK. While a short-lived 
excited state is desirable for quickly slowing a molecular beam, in order to decrease 
the final temperature of the MOT a narrower linewidth transition is needed [25] [22]. 
After careful study of the level structure of YO for the primary cooling transition, 
we identify a transition in YO that will enable narrow line cooling to the Doppler 
temperature Tp = H r )/2kB ~ 5 pK where h is the reduced Planck constant, 7 is 
the linewidth of the transition, and kp Boltzmann’s constant. Our proposed scheme 
will only require the addition of one laser to the primary cooling setup of the YO 
molecule. Narrow line cooling will greatly decrease the temperature limit of direct 
molecular cooling processes, and allow the loading of molecules into an optical dipole 
trap for possible evaporative cooling [30] . In this paper, we estimate the lifetime of 
the A' 2 A 3/2 state and calculate its vibrational branching rates in order to characterize 
the transition and assess its feasibility for use in a MOT. Additionally, we detail the 
cycling transition that will be used to create a narrow-line MOT and simulate its 
effects and cooling trajectory. 

2. Cooling process 

The entire process to cool YO from initial production down to the pK scale involves 
a number of steps, as depicted in Figure CD We form YO via laser ablation of 
a Y 2 O 3 target. The initial temperature after ablation is high both rotationally 



Figure 1: Cooling process schematic. After production via ablation (a), YO molecules 
thermalize with helium at 3.5 K in a two stage buffer gas cell (b). Upon exiting, 
the beam has velocity ~70 m/s. Longitudinal white light slowing produces molecules 
<10 m/s (c), loadable into a MOT. The orange MOT will produce cold molecules less 
than 1 mK (d). Narrow line cooling can then be performed on the red transition (e). 

and translationally, whereas any optical cycling scheme requires a large number of 
molecules in a single state. To meet this requirement, we utilize a two-stage cryogenic 
buffer gas cell m to produce cold YO molecules. The buffer gas cell is anchored 
to a 3.5 K cold plate and helium gas is constantly flowed in to maintain a high 
density (~ 10 16 /cm -3 ) of helium in the cell. The helium collisions with the cell walls 
cool the helium, and YO collisions with the helium buffer gas then rotationally and 
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translationally cool the YO molecules. The second stage of the cell has sufficiently 
low helium density to ameliorate boosting of the YO beam due to collisions with the 
helium gas at the aperture. The net result is a beam of YO molecules with forward 
velocity 70 m/s with a 40 m/s full width at half maximum, corresponding to a ~7 K 
longitudinal temperature. Longitudinal slowing of the YO beam is accomplished |21j 
with a combination of frequency chirping and white light modulation of the counter- 
propagating slowing laser, resulting in a population of YO molecules with speeds 
under 10 m/s, which are sufficiently slow to load into a first stage (orange) MOT. 
Further cooling can be accomplished with a narrow-line MOT, analagous to methods 
demonstrated for several atomic species [291 [321HH [3U [351HH |37] . 

3. Level structure of YO 

The molecule YO has three electronic states that are of primary concern, the X 2 E 
ground state, the A 2 A 1 / 2 excited state, and the A' 2 A 3 / 2 intermediate state, as shown 
in Figure [U(a). For describing the angular momentum coupling in these states we 
follow the notation used by Brown and Carrington [38] ■ When discussing a specific 
value of angular momentum, the factor of h is omitted. For the states we consider in 
YO, we have electronic spin 5 = 1/2 and nuclear spin 1= 1/2. Electron orbital angular 
momentum is labeled L, while rotational and total angular momenta are labeled R 
and F respectively. Additionally, we use intermediate angular momenta to describe 
the states: J = F — I, and N = J S = R + L. Parity of states is indicated by (+) 
or (-). As well as possessing an array of rotational levels, each state also has a ladder 
of vibrational levels labeled with quantum number v. 

We begin by considering the electronic ground state X 2 E, which can be described 
by Hund’s case (b^s). The effective Flamiltonian that describes the coupling is: 

H eS = 7 S • N + 61 • S + cI z S z 

where 7 = —9.2 MHZ, b = —763 MHz, and c = —28 MHz, representing the spin- 
rotation, Fermi contact, and electron-nuclear spin dipolar interaction respectively [39]. 
The strong Fermi contact interaction couples the electron and nuclear spin to form 
intermediate angular momentum G = I + S. The total angular momentum is then 
formed by the additional coupling of N through the spin-rotation interaction yielding 
F = G + N. We also note for E states, orbital angular momentum L is zero and so 
N = R. For A>0, levels can be expressed as a superposition of J levels J=N ±1/2 
(J= 1/2 for N= 0). 

The excited A 2 n 3 / 2 and A , 2 A 3 / 2 states can both be described by Hund’s case (a). 
The rotational states are labeled by angular momenta J and F, though the hyperfine 
interaction in these states is small and is not resolvable via optical transitions with the 
Y 2 E. The quantum number 12 is given by the projection of L + S onto the internuclear 
axis. For n states, L= 1, while for A states L= 2. So, for the A 2 n states, 12 is 1/2 or 
3/2, and for A' 2 A 12 is 3/2 or 5/2. The value of 12 is indicated by the subscript on 
the molecular term symbol, e.g. A' 2 A 3 / 2 has 12 =3/2. Additionally, the two senses of 
projections of L onto the internuclear axis give rise to a pair of states with opposite 
parity for each J level. For the A 2 A 1 / 2 state, the lambda doubling term p is 4.5072 
GHz, indicating a splitting of 4.5 GHz for the opposite parity levels in the J = 1/2 
level [40j. In the A' 2 A 3 / 2 state, the lambda doubling is expected to be small, and has 
so far not been measured. 
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The X 2 Yi to transition (u”=0 to v'=0) is at 614 nm, as depicted in 

Figure[2|a). While ^4 2 n 1 / 2 primarily decays back to A' 2 £, it has a ~3xl0~ 4 chance 
of decaying to A' 2 A 3 / 2 [II]- As will be discussed in more detail later, while the A! 2 A 3 / 2 
to A 2 £ transition is electric dipole forbidden, it does in fact occur due to mixing of 
the the nearby A 2 n 3 / 2 state with the A' 2 A 3 / 2 state [41j . This narrow transition has 
a wavelength of 690 nm. 
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Figure 2: Relevant level structure of the YO molecule, (a) Electronic structure and 
ground state rotational structure. Solid lines indicate the broad (orange) and narrow 
(maroon) transitions that will be used to make MOTs in YO. (b) Vibrational structure 
and repumping scheme. Solid lines indicate lasers for first stage cooling (orange) and 
vibrational repumps (red). Dashed lines show spontaneous decays and are labeled 
with Franck-Condon factors (c) Ground state rotational and hyperfine structure. For 
each G level, F ranges from |G — N\ to |G + iV|. 


4. First stage laser cooling 

Our first stage of cooling has a main cycling transition at 614 nm from A 2 £, N”=1 
to A 2 n^/2, J= 1/2, as shown in Figure [2](a). Optical cycling has been demonstrated 
on this transition as a one and two dimensional MOT [221 . and for molecular beam 
slowing [2T]- Decays to rotational levels other than IV” =1 are forbidden by parity and 
angular momentum selection rules. We address vibrational branching by repumping 
the v"=l level directly to the x4 2 n 1 / 2 ^ w ,= 0 level using 648 nm light (Figure IU(b)). 
The v"—1 level is repumped using a two step process, first exciting it to A 2 IIj/2, 
v'=l at 649 nm, which predominantly decays to the v” = 1 level, and then allowing 
the v"=l repump light to take the molecule back to the v”=0 level. Molecules that 
decay to X 2 £ via A' 2 A 3 / 2 end in the 7V”=0 or 2 levels because of parity and angular 
momentum rules. To maintain cycling, we apply microwave mixing resonant with the 
iV”=0y->-iV”=l and the !V”=leiV”=2 rotational transitions, depicted as gray arrows 
in Figure [2jc). In total, we estimate that the setup yields a closed cycling transition 
up to the ~ 10 -6 level of enclosure, sufficient to slow longitudinally and trap a subset 
of the molecules. 

In addition to being able to perform Doppler cooling on the X 2 Y, to A 2 I \ l / 2 
transition, the G”=l: F” =1,2 levels are subject to Zeeman splitting in an applied 
magnetic field, allowing the use of a MOT. For this transition, the greater multiplicity 
of states in the ground state means that coherent dark states can form |42j . limiting 
the total scattering rate. In our setup, we avoid this issue by switching our laser 
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Table 1: Identification of A' 2 E to A' 2 A 3 / 2 transitions 


Transition 

Frequency (THz) 

N”=0,G”=1 ->■ A ,2 A 3/2 (-) 
N”=1.G”=1 -y A ,2 A 3/2 (+) 
N”=2,G”=1 -y A ,2 A 3/2 (-) 

434.70310(3) 

434.67989(3) 

434.63333(3) 


polarization and the magnetic field synchronously at a rate approximately equal to 
the transition linewidth, eliminating long-lasting dark states. Type-II MOTs like this 
one (operating from F" to F' = F" — 1) tend to be weaker and larger than their type-I 
counterparts (operating from F” to F' = F” + 1), with sizes on the order of a mm 
[43] . We anticipate the final temperature of the orange MOT to be on the order of 
1 mK. This is larger than the Doppler temperature Tjj = 110/iff, but typically type-II 
MOTs do not reach the Doppler limit [44], [23]. However, the rapid polarization and 
magnetic field switching should negate the weakness of the type-II MOT in our case 
[45] and may allow a lower final temperature of the orange MOT. 

5 . Spectroscopy of the A' 2 A 3 / 2 state 

In order to characterize the properties of a second stage red MOT on the AT 2 E to 
A! 2 A 3/2 transition, we first need an estimate of the linewidth of the transition. We 
acquire this estimate by comparing the absorption strength of the red transition to that 
of the known orange transition. To that end, we have identified the IV”=0,1,2:G”=1 
to A! 2 A 3/2 transitions using absorption spectroscopy in our cryogenic buffer gas cell, 
as shown in Table [T] in agreement with previous work |46] . 

The energy splitting between N”=0—> A’ 2 A 3 / 2 and IV” =2 —> A’ 2 A 3 / 2 is IV”(IV”-|- 
1)B 0 , with N” = 2, where Bq is the rotational constant for A 2 E, t>”=0. Evaluating 
this yields Bq = 11.61 GHz, which agrees with the known rotational splitting [41] . 
By comparison of the absorption on the IV” =1 —> A , 2 A 3 / 2 transition with the IV” =1 
-s-4 2 n 1/2 transition and the known 33 ns lifetime of the 4 2 n 1 / 2 state[37], we estimate 
that the lifetime of the A' 2 A 3 / 2 state is r =425 ns. We note that the IV” = 0,2 
to A! 2 A 3 / 2 transitions address a different parity level than the IV” =1 to A' 2 A 3 / 2 
transition, but we do not find a different value of H, indicating that the parity splitting 
of the 4' 2 A 3 / 2 state is smaller than the Doppler broadened linewidth of ~ 60 MHz 
inside the cell. 

6. Mixing of the A ,2 A 3 / 2 and A 2 n states 

In addition to performing an absorption measurement, we also calculate a lifetime 
estimate. While electric dipole forbidden, the decay from the A , 2 A 3 / 2 state is almost 
entirely due to mixing with the A 2 n 3 / 2 state, on the order of 3.5% [IT]. The 
molecular potential can be approximated using the Rydberg-Klein-Rees (RKR) [48] 
method for the states of interest. The curves for the A 2 E and A 2 H states can 
be corroborated with previous work [49]. After using a finite difference method to 
determine the wavefunctions to evaluate the vibrational overlap integrals and the 
{v\l/r 2 \v') integrals we use the procedure outlined by Chalek and Gole [IT] to 
compose the Hamiltonian including the electronic, rotational, vibrational, spin orbit, 
and coriolis operators. By diagonalizing this Hamiltonian for J= 3/2, we express the 
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Table 2: Vibrational branching rates for various transitions of interest. 


State 1 

State 2 

Vl 

v 2 

Qv^v 2 

State 1 

State 2 

Vl 

v 2 

Qviv 2 

^4 2 IIi/2 

X 2 S 

0 

0 

0.9924 

7 l 2 IIl /2 

A 2 £ 

1 

0 

0.0076 

a 2 u 1/2 

A' 2 £ 

0 

1 

0.0074 

7F1I1/2 

A 2 £ 

1 

1 

0.9752 

xt 2 n 1/2 

X 2 S 

0 

2 

2.4xl0~ 4 

a 2 yi 1/2 

A 2 £ 

1 

2 

0.0164 

a 2 u 1/2 

A' 2 £ 

0 

3 

2xl0 -6 

^4 2 IIi/2 

A 2 £ 

1 

3 

7.8x 10 -4 

a 2 u 3/2 

^ 4 /2 A3/2 

0 

0 

0.9078 

A ' 2 A 3/2 

A 2 £ 

0 

0 

0.9422 

^ 2 n 3 / 2 

^ 4 /2 A3/2 

0 

1 

0.0860 

A ' 2 A 3/2 

A' 2 £ 

0 

1 

0.0568 

^ 2 n 3/2 

A ' 2 A3/2 

0 

2 

0.0059 

A ' 2 A 3/2 

A 2 £ 

0 

2 

0.0011 

- 4 2 n 3 / 2 

A 2 A 3/ 2 

0 

3 

3.3x 10~ 4 

A' 2 A 3/2 

A 2 £ 

0 

3 

<10~ 4 


A' 2 A 3 / 2 eigenstate in terms of the Born-Oppenheimer states, the first few terms of 
which are: 


tfA3 /a ,»'=o> « 0.9865 |A 3/2 y = 0) - 0.1582 \n 3/2 ,v' = 0) 

-0.0369 |n 3/2 , v' = 1) + 0.0176 |n 1/2 , v' = 0) + ... 


These coefficients squared indicate the mixing of the A' 2 A 3 / 2 ,i/=0 state with 
other states, with the largest being a 2.5% mixing with the A 2 n 3 / 2 , v'= 0 state. 
Summing over the various Born-Oppenheimer states and their associated calculated 
Franck-Condon factors yields estimates for the vibrational branching ratios for the 
A! 2 A 3 j 2 to V 2 E transition. These and other state branching ratios are compiled in 
Table [2] 

Using 2.5% mixing and the known ~30 ns lifetime of the A 2 n 3 / 2 state [47], we 
estimate the lifetime of the A' 2 A 3 / 2 state as 1200 ns. These two estimates of the 
lifetime of the A' 2 A 3 / 2 state, 425 ns from the absorption measurement and 1200 
ns from calculation are of similar order and give us a starting point for further 
calculations. From here on, we use an intermediate estimate of r= 1 /xs, yielding 
a linewidth of 7 = 2-7TX160 kHz. 

7. Narrow-line cooling on the V 2 E to A' 2 A 3 / 2 transition 

Once molecules have been cooled with the orange MOT, we wish to further cool them 
utilizing the narrow V 2 E to A' 2 A 3 / 2 transition. This is acheivable by exciting from 
the V” = 0 level to the A' 2 A 3 / 2 , J'=3/2 level, as in Figure^a). Decay from A' 2 A 3 / 2 
ends in either TV” =0 or 2. To maintain cycling, we use microwave mixing to address the 
i>”=0, TV”=0-f->T, and TV”=le->-2 rotational transitions, with frequencies 2B$ and ABq. 
Similarly, higher vibrational levels can be mixed by applying microwaves resonant 
with the TV” =0-0-1 and TV”= 1 - 0-2 rotational transitions of that vibrational level, e.g. 
at 2B\ and AB\ for u”=l. Additionally, this microwave mixing scheme allows fo’>0 
vibrational levels to be repumped with the same lasers used in the orange MOT. 
The rotational constants for the first few vibrational levels have been determined and 
so mixing will be straightforward to implement m isi]- The saturation intensity, 
I sat — (7t/ic)/(3A 3 t), is 60 /iW/crn 2 , which is obtainable with commercially available 
diodes. 

In addition to the possibility of an optical molasses based on this transition, it is 
apparent that a narrow-line MOT is also achievable. A MOT can be made using any or 
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Figure 3: Level structure for the narrow-line MOT. (a) Maroon arrow indicates the 
transition used in the red MOT. Gray lines indicate microwave mixing which closes 
rotational branching. Red line indicates r>”=1 repump, the same as used in the 
orange MOT. (b) Zeeman structure and MOT scheme, g represents the corresponding 
magnetic Lande factor. Large arrows indicate counter propagating circularly polarized 
laser beams. Both F”= 0 —> F '=1 and F”=l— >F'=2 transitions can effect a MOT 
action (maroon arrows). 


all F” states due to a non-zero g-factor for both the ground and excited states, but the 
number of frequency components and polarization considerations go up considerably 
with the number of states {45j. The simplest possible MOT action in this system is 
for the transition from the iV”=0 level to the A ,2 A 3 /2 state. The IV” =0 level has G 
values 0 and 1, split by ~760 MHz. A type-I MOT can be performed on both sublevels, 
from G”=F”=1 to F'= 2 and from G” = F” = 0 to F'= 1, as depicted in Figure E^b). 
Additionally, since in type-I MOTs no Zeeman dark states can form, experimental 
implementation is simplified because it will not be necessary to perform the magnetic 
field switching and polarization chopping used in the orange MOT. 

8. Narrow-line MOT simulation 

We simulate the red MOT in one dimension with counter-propagating circularly 
polarized light with components addressing both the IV” =0, G”=0 and G”=l 
manifolds and exciting to available A r2 A 3 /2 states. Initial conditions are given by 
a 1 mK distribution using the parameters of an orange MOT with a trap oscillation 
frequency of ujmot = 27rxl55 Hz [22]. Because the Doppler broadened linewidth 
at 1 mK is significantly greater than the narrow transition, we model the laser light 
with multiple frequency components spanning the broadened linewidth, each with 
saturation parameter so [29], .32 . We use a carrier plus two pairs of sidebands spaced 
by 27rxl MHz, with the detuning of the highest energy sideband from resonance 
labeled <5i. Since microwave mixing and vibrational repumping are both fast compared 
to the optical cycling on this transition, we make the following assumptions about 
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the population distribution in the ground state. Because of the 760 MHz hyperfine 
splitting, G”=0 manifolds are only mixed with other G”=0 states, and G”=l manifolds 
with other G”=l states. So, population in a given G” manifold is weighted by the 
branching ratio from the A' 2 A 3/2 state, which is 1/4 and 3/4 into the G”=0 and 1 
levels respectively. For a given G” level, the population is evenly distributed among 
all available u”=0, N” =0,1,2 states due to the microwave mixing. For a two level 
system, the force from a laser beam frequency component resulting in absorption and 
followed by spontaneous emission can be written as [521 

_ hkso"f/2 

~ 1 + so + (2<V7) 2 

where k is the wavenumber, So is the saturation parameter, S the detuning experienced 
by the molecule, and 7 = 1/r. In a MOT, S can be written as 5 = 60 — kv + p, Ar/h, 
where <5o is the nominal detuning of the laser frequency component from resonance, 
v the velocity, A the magnetic field gradient, r the molecule’s position, and /.1 the 
differential magnetic moment for the transition. To compute the average force felt by 
an TV” = 0 molecule, the force is averaged over the transitions available to all N” = 0 
sublevels, and summed over both laser beams and all frequency components. Since 
we only optically address iV”=0 which contains 1/9 of the microwave mixed states, 
the average force felt by a molecule mixed between N” =0,1,2 is Ftotai=FN”= 0 / 9 - 
In addition, we simulate heating due to spontaneous emission by adding a random 
momentum kick when a photon is absorbed. For small detunings, the force reduces 
to F tota i = —/3v — k r, where /3 is the damping cofficient and k the spring constant, 
which together characterize the trap. 


(a) „ (b) (c) (d) 



Figure 4: Results of simulation, (a) Velocity dependent force due to optical molasses. 

(b) Position dependent force due to magnetic field gradient for stationary molecules. 

(c) 1 mK distribution produced from orange MOT (blue) and distribution after red 
MOT has been on for 20 ms (red), (d) Temperature vs time after red MOT is turned 
on and the orange MOT is turned off. 

We first model an optical molasses as shown in Figure QJa) for the parameters 
<5i = -I .57 and so = 3, chosen to provide force to all moving molecules in an initial 
1 mK distribution. From the slope around zero velocity, we can extract /? and so 
the damping rate Td = 2/3/m = 950 Hz. Likewise, we extract n from the rate of 
change of force with position for zero velocity molecules, as shown in Figure SKb). 
Under the same parameters, and magnetic field gradient A = 3 G/cm, we find a MOT 
oscillation frequency u>mot = \J «/ m °f 27 tx33 Hz. The characteristic trapping time is 
then Te,/uj\ iot = 20 ms [52!. Figures|U(c) and|U(d) show the results of a Monte-Carlo 
simulation of the narrow-line MOT. By comparing the initial temperature distribution 
(blue) and after 20 ms (red), we see clear phase space compression. The temperature 
after 20 ms is less than 10 /j,K . The simulation indicates a photon scatter rate ~10 4 
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per second, and cooling to <10 pK requires only ~200 photons. A naive worst case 
extrapolation to three dimensions is yielded by dividing n and /? by three, giving 
uj 3 e> = 27tx19 Hz, r 3£ > =317 Hz, and requires ~600 photons to cool. We estimate the 
v v = 3 limited MOT lifetime to be >1 second. We note that for a molecule at rest, the 
maximum acceleration due to the three dimensional MOT is ^5 times that of gravity, 
allowing successful trapping. In sum, the simulation indicates that the red MOT is 
capable of trapping and further cooling the molecules captured in an orange MOT, 
requiring only the addition of one 690 nm laser with frequency components split by 
760 MHz, obtainable with an acousto-optic modulator, to the first stage system, one 
for the G”=0 and the other for the G"=1 to A' 2 A 3 / 2 transition. 

9. Conclusion 

We have identified a transition in YO suitable for secondarily cooling a MOT on the 
X 2 Y, to A 2 n 1/2 transition by using a narrow-line transition from the A 2 E to A ,2 A 3 / 2 
state. Two estimates indicate a lifetime of r «1 /is for the A' 2 A state, allowing a final 
temperature of ~10 pK. We have also determined that the vibrational branching from 
the A' 2 A 3 / 2 state is favorable enough to not necessitate any more vibrational repump 
lasers than the first stage MOT. Simulations indicate parameters for the red MOT of 
a damping rate of 317 Hz and a MOT frequency of 27tx19 Hz. 
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